The Kondo Box: A Magnetic Impurity in an Ultrasmall Metallic Grain 
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We study the Kondo effect generated by a single magnetic impurity embedded in an ultrasmall 
metallic grain, to be called a "Kondo box" . We find that the Kondo resonance is strongly affected 
when the mean level spacing in the grain becomes larger than the Kondo temperature, in a way 
that depends on the parity of the number of electrons on the grain. We show that the single- 
electron tunneling conductance through such a grain features Kondo-induced Fano-type resonances 
of measurable size, with an anomalous dependence on temperature and level spacing. 

PACS numbers: 72.15.Qm, 73.20.Dx, 73.23.Hk, 71.27. +a 



What happens to the Kondo effect when a metal sam- 
ple containing magnetic impurities is made so small that 
its conduction electron spectrum becomes discrete with a 
non-zero mean level spacing A? More specifically, when 
will the Kondo resonance at the Fermi energy ep that 
characterizes bulk Kondo physics begin to be affected? 
This will occur on a scale A ~ Tjc, the bulk Kondo tem- 
perature, since a fully-developed resonance requires a fi- 
nite density of states (DOS) near ep, and A will act as 
low-energy cut-off for the spin scattering amplitude. 

To achieve A > Tk, the sample would have to be an 
ultrasmall metallic grain containing magnetic impurities, 
to be called a "Kondo box": for example, for a metal- 
lic grain of volume (15nm) 3 -(3nm) 3 and kp ~ lA , 
the free-electron estimate A = 1/JVo — 2-k 2 K 2 / 1 (mkpVol), 
with N the bulk DOS near e F , gives A ~ 0.5 - 60K, 
which sweeps a range including many typical Kondo tem- 
peratures. The discrete DOS of an individual grain of 
this size can be measured directly using single-electron 
tunneling (SET) spectroscopy Jl|,||, as shown by Ralph, 
Black and Tinkham Q in their studies of how a large 
level spacing affects superconductivity. Analogous ex- 
periments on a Kondo box should be able to probe how 
a large A(~ Tk) affects Kondo physics. 

In this Letter we study this question theoretically. We 
find (1) that the Kondo resonance splits up into a series 
of sub-peaks corresponding to the discrete box levels; (2) 
that its signature in the SET conductance through the 
grain consists of Fano-like line shapes with an anomalous 
temperature dependence, estimated to be of measurable 
size; (3) an even/odd effect: if the total number of elec- 
trons on the grain (i.e. delocalized conduction electrons 
plus one localized impurity electron) is odd, the weight 
of the Kondo resonance decreases more strongly with in- 
creasing A and T than if it is even. 

The model: — For the impurity concentrations of 0.01% 
to 0.001% that yield a detectable Kondo effect in bulk 
alloys, an ultrasmall grain of typically 10 4 — 10 5 atoms 
will contain only a single impurity, so that inter-impurity 
interactions need not be considered. We thus begin by 



studying the local dynamics of a single impurity in an 
isolated Kondo box, for which we adopt the (infinite U) 
Anderson model with a discrete conduction spectrum, in 
the slave-boson representation: 

H = H +e d J2fti* + vJ2( c ]<r bi f° + h - c -)' (!) 

where Hq = ^22ija £ i^ja c ja- Here a denotes spin and 
the cl a create conduction electrons in the discrete, de- 
localized eigenstates \ja) of the "free" system (i.e. with- 
out impurity). Their energies, measured relative to the 
chemical potential /i, are taken uniformly spaced for sim- 
plicity: £j = jA + Eq — fj,. As in S, we follow the 
so-called orthodox model and assume that the e,-'s in- 
clude all effects of Coulomb interactions involving delo- 
calized electrons, up to an overall constant, the charg- 
ing energy Eq- The localized level of the magnetic 
impurity has bare energy Sd far below ep, and is rep- 
resented in terms of auxiliary fermion and boson op- 
erators as d\ — f^b, supplemented by the constraint 
So- Safe + b^b = 1 [Q, which implements the limit 
U — * oo for the Coulomb repulsion U between two elec- 
trons on the d-leve\. Its hybridization matrix element 
v with the conduction band is an overlap integral be- 
tween a localized and a delocalized wave-function, and, 
due to the normalization of the latter, scales as ( Vol)^ 1 ! 2 . 
Thus the effective width of the ci-level, T = irv 2 / A, is 
volume-independent, as is the bulk Kondo temperature, 
Tk = \j2YDjix exp(— ned/^T), where D is a high energy 
band cutoff. To distinguish, within the grand canonical 
formalism, grains for which the total number of electrons 
is even or odd, we choose fj, either on (/i = £q) or half- 
way between two (fj, = Eq + A/2) single-particle levels, 
respectively ||. 

NCA approach: — We calculated the spectral den- 
sity Adaito) of the impurity Green's function Gd&(t) = 
— i8(t){{d a (t), dj.(0)}) using the noncrossing approxima- 
tion (NCA) ||. For a continuous conduction band, the 
NCA is known to be reliable down to energies of O.ITk 
or less, producing spurious singularities only for T below 
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FIG. 1. Impurity spectral function Ada(uj) for various 
values of A at T — 0.5Tk; (a) even, (b) odd total number of 
electrons. The individual curves are offset by one unit each. 

this scale HQ . Since these are cut off by the level spacing 
A in the present case, we expect the NCA to be semi- 
quantitatively accurate over the entire parameter range 
studied here (T and A between 0.1 and 5Tk)- Denoting 
the retarded auxiliary fermion and boson propagators by 
G ftF {u) = (uj-s d - E/aH)- 1 , G b (u) = {u- E^))" 1 , 
respectively, the selfconsistent NCA equations read 

E /CT H=T J ^.[l-f(e)]A^(e)G b (cj-e), (2) 

E 6 (uO / ^mA<V{e)G f „(u, + E), (3) 

where f(u>) — l/[exp(w/T) + 1]. The finite grain size 
enters through the discreteness of the (dimensionless) 
single-particle spectral density of the box without im- 
purity, A^J(uj) = AJ2jS{uJ — ej). (We checked that 
all our results are essentially unchanged if the Dirac 
8's are slightly broadened by a width 7 < 0.12V.) In 
terms of the auxiliary particle spectral functions Af^ = 
— ilmG/,b, Ada(u) is given by (for details see 0) 

A da (u>) = f de[c-^ + c-^-^]A fa (e)A b (e - u) . (4) 



Numerical results: — The results obtained for A da (cj) 
by numerically solving the NCA equations (0) to (0) 
for various T and A are summarized in Figs, [j] and|2. 
(We have checked that the equation-of-motion method 
|| yields qualitatively similar results for all quanti- 
ties discussed below.) For A -c T, the shape of the 
Kondo resonance is indistinguishable from the bulk case 
(A — » 0); when A is increased well beyond T, how- 
ever, it splits up into a set of individual sub-peaks. 
With decreasing temperature (at fixed A), each sub- 
peak becomes higher and narrower; its width was found 
to decrease without saturation down to the lowest tem- 
peratures for which our numerics were stable (T ~ 
0.2 A). This agrees with the expectation following 
from the Lehmann representation at T = 0, A da {uj) 
= £„ [|(n|4|0)| 2 5(a; - A„) + \(n\d a \0)\ 2 6(cj + A n ) , 
namely that the sub-peaks should reduce to <5-functions 
with zero width, located at the exact excitation energies 
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FIG. 2. Even/odd dependence of the average peak height 
of the Kondo resonance, as function of T/Tk- For an even 
box (o), we averaged Ada over a range A centered on its cen- 
tral sub-peak, for an odd box (•) over a range 2A centered on 
its central two sub-peaks (as indicated by arrows in Fig. 1). 
The inset shows the same quantity as function of A/Tk- Nu- 
merical uncertainties are smaller than the symbol sizes. 

A„ = E n — Eo of the full Hamiltonian (whose spectrum 
will be discrete too, with mean level spacing of the same 
order as A, as follows from the exact finite-size solutions 
of the Kondo model || ) . 

Despite developing sub-peaks, the Kondo resonance re- 
tains its main feature, namely significant spectral weight 
within a few Tk around the Fermi energy, up to the 
largest ratios of A/max(T, Tk) (— 5) we considered. This 
implies that the Kondo correlations induced by the spin- 
flip transistions between the (i-level and the lowest-lying 
unoccupied j-levels persist up to remarkably large values 
of A/max(T, Tjc) j[o|. However, they do weaken system- 
atically with increasing A, as can be seen in the inset 
of Fig. 2, which shows the average peak height of the 
Kondo resonance (which quantifies the "strength" of the 
Kondo correlations) as function of A at fixed T: the peak 
height drops logarithmically with increasing A once A 
becomes larger than about T. Conversely, at fixed A, it 
drops logarithmically with increasing T once T becomes 
larger than about 0.5 A (main part of Fig. 2), thus repro- 
ducing the familiar bulk behavior. Qualitatively, these 
features are readily understood in perturbation theory, 
where the logarithmic divergence of the spin flip ampli- 
tude, t(u>) oc Y" 1 . , is cut off by either T or A, 
whichever is largest. 

Parity Effects: — For A ^> T, the even and odd spec- 
tral functions A da in Fig. 1 differ strikingly: the former 
has a single central main peak, whereas the latter has two 
main peaks of roughly equal weight. This can be under- 
stood as follows: For an even grain, spin-flip transitions 
lower the energy by roughly Tk by binding the d elec- 
tron and the conduction electrons into a Kondo singlet, 
in which the topmost, singly-occupied j level of the free 
Fermi sea carries the dominant weight, hence the single 
dominant peak in Ada-- For an odd grain, in contrast, the 
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free Fermi sea's topmost j level is doubly occupied, block- 
ing such energy-lowering spin- flip transitions. To allow 
the latter to occur, these topmost two electrons are redis- 
tributed with roughly equal weights between this and the 
next-higher-lying level, causing two main peaks in Ada 
and reducing the net energy gain from Tk by an amount 
of order A. This energy penalty intrinsically weakens 
Kondo correlations in odd relative to even grains; indeed, 
the average Ada peak heights in Fig. 2 are systematically 
lower in odd than in even grains, and more so the larger 
A and the smaller T. 

SET conductance: — The above physics should show 
up in SET spectroscopy experiments: When an ultra- 
small grain is connected via tunnel junctions to left (L) 
and right (R) leads |llf and if the tunneling current 
through the grain is sufficiently small (so that it only 
probes but does not disturb the physics on the grain), 
the tunneling conductance GiV) as function of the trans- 
port voltage V has been demonstrated jjj to reflect the 
grain's discrete, equilibrium conduction electron DOS. 
Such measurements are parity-sensitive [0 even though a 
non-zero current requires parity-fluctuations, since these 
can be minimized by exploiting the huge charging ener- 
gies (E c > 50K) of the ultrasmall grain. To calculate the 
SET current, we describe tunneling between grain and 



leads by H t = Y,kia a ( u k 
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creates a spin a electron in channel k of lead a £ {L, R}. 
Neglecting non-equilibrium effects in the grain, the tun- 
neling current has the Landauer-Buttiker form |l2| 



duj Fy (w) 
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Actios), (5) 



where F v {u) = f{u - eV/2) 
spectral density of G C j a {i) 
and 7 ?=27r£ fc 



l kja I 



- f(w + eV/2), A c . ja is the 
= -tf(tX{c^(t)/ C t a (0)}), 
Neglecting the aja depen- 
dence of 7, the current thus is governed by the conduc- 
tion electron DOS, A c (ui) — Ylja-A c j a (u)). Exploiting a 
Dyson equation for G c j a , it has the form 

A c {lo) = --^Im [GojctM + v 2 [G%{u>)] 2 G d a{u>) 



3 a 



where G^L = 1/(uj—e 
tron Green's function 



,o_riO + ) is the free conduction elec- 
[13] , and the corresponding Kondo 
contribution to the conductance G(V) — dI{V)/dV — 
G (V)+6G(V) is 



e 2 r 

5G{V) = -W- 
n 7r 



UJ) X 



(0) 



F v {u) 



F v (e 3 ) dF v (uj)/dLo 



UJ — Ej 



withFy(w) = -~-£(f(uj-eV/2) + f(uj + eV/2))/2. Even 
though Kondo physics appears only in the subleading 
contributions to A c (lu) and G(V), these are proportional 
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FIG. 3. Differential conductance of (a) an even and (b) 
and odd Kondo box for A = 4Tk- Dotted, dashed and 
solid lines give, respectively, the bare conductance Gn(V), 
the Kondo contribution SG(V), and the total conductance 
Go + 5G. The curves with larger or smaller amplitudes corre- 
spond to T = 0.7Tk and ITk, respectively; all are normalized 
such that the average of the bare conductance Go(V) = 1. 



to v 2 = rA/7r and thus grow with decreasing grain size. 

SG(V) and G(V) are shown in Fig. || and have rather 
irregular structures and line-shapes. The reason for this 
lies in the interference between Gda and [G^j^] 2 in A c , 
and correspondingly between Ada and the bracketed fac- 
tor in (H) for 8G(V). This interference is reminiscent of 
a Fano resonance Jl5|], which likewise arises from the in- 
terference between a resonance and the conduction elec- 
tron DOS. Incidentally, Fano-like interference has been 
observed in STM spectroscopy of a single Kondo ion on 
a metal surface [ |l6fl , for which the conduction electron 
DOS is flat. In contrast, for an ultrasmall grain it con- 
sists of discrete peaks, reflected in the last factor in Eq. 
(^). This leads to a much more complex interference 
pattern, which does not directly mirror the specific peak 
structure of Ada (w) discussed above. 

Nevertheless, G(V) does bear observable traces of the 
Kondo effect, in that the interference pattern shows a 
distinct, anomalous T -dependence, due to that of the 
Kondo resonance. In particular, the weights Wj under 
the individual peaks of G(V) become T-dependent. (In 
contrast, the weight Wq under an individual peak of the 
bare conductance Gq(V) is T independent, since the T 
dependence of the peak shapes of Go are determined 
solely by df{ui)/du>.) This is illustrated in Fig. [|, which 
shows the T dependence of the weights W\ and Wi of 
the first and second conductance peaks (counted relative 
to V = and labelled 1,2 in Fig. When T decreases 
at fixed A = Tk, both W\ and W% decrease, while at 
fixed A = 3Tk, W\ decreases whereas W2 increases. The 
fact that the weights can either increase or decrease with 
decreasing T results from the constructive or destructive 
Fano-like interference effects discussed above. Moreover, 
at the larger value for A, both W\ and W2 develop a 
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FIG. 4. Anomalous temperature dependence of the weights 
W\ (circles) and W2 (diamonds) of the first two conductance 
peaks of an even grain (o,0) and an odd grain (•,♦), for 
A = 3Tk and A = Tk (inset). We calculated the weights 
under G(V) over a fixed range V = A between successive 
minima of the lowest-T 5G(V) curve (see arrows in Fig. 3). 

parity effect in the strength of their T dependence. Since 
the peak weights in Fig. |] change by up to ~ 10% as 
the grain is cooled below Tk , it should be possible to ex- 
perimentally |y detect their Kondo-induced anomalous 
T-dependence. 

Coherence length: — The condition A > Tk implies a 
relation between sample volume and the much-discussed 
spin coherence length £k = 2t:vf /Tk, namely (in 3D) 
Vol < n^K/kp. Note that this relation involves both the 
small length scale 1/kF and the sample's volume, and not 
the smallest of its linear dimensions, say L. This implies 
that the length scale below which purely finite-size in- 
duced modifications of Kondo physics can be expected is 
not set by £k alone Jl7[ , and indeed may be considerably 
smaller than £k ■ This is why such modifications were not 
found in the numerous recent experiments having L < £k 
for one or two sample dimensions fll8||T9| ]. 

In conclusion, we have analysed the Kondo effect in 
an ultrasmall metallic grain containing a single magnetic 
ion. The presence of a new energy scale in the system, 
the mean level spacing A, leads to rich physical beha- 
vior when A ~ Tk, including a distinct even/odd effect. 
Our NCA calculations, which give a semi-quantitatively 
reliable estimate of the size of the effects to be expected 
in future experiments, predict that the SET conductance 
through such a grain has a Kondo-induced anomalous T- 
dependendence of up to 10%. Since the effects discussed 
in this work result, above all, from the discrete density 
of states near sf, they should be generic for ultrasmall 
grains, i.e. robust against including randomness in the 
model, like j-dependent level spacings Aj and hybridiza- 
tion matrix elements Vj . An interesting extension of this 
work would be to study the effects of a magnetic field, 
which couples not only to the magnetic impurity but also 
Zeeman-splits each discrete j-level on the grain. 
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